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This paper reports measurements of atmospheric neutrino and antineutrino interactions in the 
MINOS Far Detector, based on 2553 live-days (37.9 kton-years) of data. A total of 2072 candidate 
events are observed. These are separated into 905 contained- vertex muons and 466 neutrino-induced 
rock-muons, both produced by charged-current !/ M and 77 M interactions, and 701 contained-vertex 
showers, composed mainly of charged-current u e and V e interactions and neutral-current interactions. 
The curvature of muon tracks in the magnetic field of the MINOS Far Detector is used to select 
separate samples of z/ M and 77 M events. The observed ratio of to events is compared with the 



Monte Carlo simulation, giving a double ratio of Rffi/Rffi = 1.03±0.08 (stat.) ±0.08 (syst.). The 
Vfj. and Vfj, data are separated into bins of L/E resolution, based on the reconstructed energy and 
direction of each event, and a maximum likelihood fit to the observed L/E distributions is used to 
determine the atmospheric neutrino oscillation parameters. This fit returns 90% confidence limits of 
|Am 2 | = (1.9±0.4) x 10~ 3 eV 2 and sin 2 2# > 0.86. The fit is extended to incorporate separate and 
F M oscillation parameters, returning 90% confidence limits of |Am 2 | — |Am 2 | = 0.6t(j'g x 10 -3 eV 2 
on the difference between the squared-mass splittings for neutrinos and antineutrinos. 

PACS numbers: 14.60. Pq 



I. INTRODUCTION 

It has now been firmly established by experiment that 
muon neutrinos produced by cosmic-ray showers in the 
atmosphere undergo oscillations. The data are well de- 
scribed by Vp — > v T neutrino oscillations, and measure- 
ments of the oscillation parameters have been made by 
Super-Kamiokande (SK) MACRO 0, Soudan 2 @ 
and MINOS (fj, [j]]. The atmospheric neutrino results are 
strongly supported by long-baseline experiments, which 
observe corresponding oscillations in accelerator beams 
of muon neutrinos. Beam neutrino measurements have 
been made by K2K @, T2K § and MINOS [IM1, 
The MINOS beam data analysis, which uses a two-flavor 
model of neutrino oscillations, returns best fit values of 
|Aro 2 | = (2.32i°:os) x 10~ 3 eV 2 and sin 2 26» = 1.00_ .oe 
for the oscillation parameters fl2| . 

The MINOS experiment has performed separate mea- 
surements of antineutrino oscillations [13l - ll5| by iden- 
tifying antineutrino interactions in the Fermilab NuMI 
accelerator beam [l6[ . A precision measurement of these 
oscillations has been made by operating the NuMI beam 
in a 77^-enhanced configuration. Using the 77^-enhanced 
data set, the v ^ oscillation parameters are measured to 
be |Am 2 | = [2.62±g:|| (stat.) ± 0.09 (syst.)] x 10" 3 eV 2 
and sin 2 20 = 0.95^;^ (stat.) ± 0.01 (syst.) [H]. Such 
studies are of interest, as an apparent difference between 
the and oscillation parameters could indicate new 
physics. In particular, separate and F M measurements 
can be used to study models of non-standard neutrino 
interactions (l7l Il8|, a nd probe CPT symmetry in the 
neutrino sector [19j, [20( . 

The SK experiment has also studied oscillations in at- 
mospheric neutrinos and antineutrinos. Although SK 
cannot distinguish from v ^ on an event-by-event ba- 
sis, they have performed a statistical analysis of their 
data and the results are consistent with equal and F„ 
oscillation parameters [2Tj . 

The MINOS experiment is able to study atmo- 
spheric neutrinos and antineutrinos separately using its 
5.4kton Far Detector, which is located 705 m under- 
ground (2070 m water-equivalent) in the Soudan mine, 
Minnesota. At this depth, the incident flux of cosmic- 
ray muons is reduced by a factor of 10 6 relative to the 
surface. By applying a series of selection requirements, 
the cosmic-ray muon background can be reduced by a fur- 
ther factor of 10 6 , yielding a clean sample of atmospheric 
neutrino signal events. The MINOS Far Detector is mag- 



netized, which enables atmospheric + N (JT ± X 
and + N — > fi + + X charged-current (CC) interactions 
to be separated based on the curvature of the muons. 

MINOS has been collecting atmospheric neutrino data 
since 2003 and has previously published charge-separated 
analyses of contained-vertex muons [6| and neutrino- 
induced rock- muons 0, [HJ, based on 418 and 854 live- 
days of data, respectively. Both samples are largely com- 
posed of atmospheric neutrino and CC interactions, 
which are identified by the presence of a primary muon 
track in the reconstructed event. For contained-vertex 
muons, the reconstructed interaction vertex is contained 
inside the fiducial volume of the detector. The sample in- 
cludes both fully-contained muons, which stop in the de- 
tector, and partially-contained muons, which exit the de- 
tector. The muon track is typically accompanied by some 
vertex shower activity, generated by the hadronic system, 
which is used to fully reconstruct the neutrino energy. 
For neutrino-induced rock-muons, the reconstructed ver- 
tex is outside the fiducial volume. The selected muons 
enter the detector in an upward-going or horizontal di- 
rection and can be either stopping or through-going. 

The atmospheric neutrino analysis presented here is 
based on an updated data set of 2553 live-days, collected 
between August 2003 and March 2011. The contained- 
vertex muon and neutrino-induced rock-muon samples 
have been combined into a single analysis, along with an 
additional sample of contained-vertex showering neutri- 
nos, which are mainly composed of v e and v e CC inter- 
actions and neutral-current (NC) interactions. The data 
are compared to the hypothesis of — > v r and v ^ — > ~D T 
two-flavor vacuum oscillations. 



II. THE MINOS FAR DETECTOR 

The MINOS Far Detector [23] is a steel-scintillator 
calorimeter, containing 486 octagonal planes of 2.54cm 
thick steel, interleaved with planes of 1 cm thick extruded 
polystyrene scintillator and air gaps of 2.4 cm thickness. 
The planes are vertical, with a height of 8 m. Each scin- 
tillator plane is divided into 192 strips of width 4.1cm, 
aligned at ±45 degrees to vertical. The direction of the 
strips alternates from plane to plane. The scintillation 
light is collected using wavelength-shifting fibers, which 
are embedded in the strips. At the ends of each strip, 
the emitted light is transported by clear optical fibers to 
multi-anode photo-multiplier tubes. 
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The detector comprises two super-modules, of length 
14.8 m and 14.0 m, separated by a gap of 1.1m. Each 
super-module is magnetized toroidally to an average field 
of 1.3 T using a current loop that runs through a 25 cm 
diameter coil hole along the central axis of the super- 
module and then returns below the super-module. The 
MINOS coordinate system is right-handed, with the y- 
axis pointed vertically upwards and the z-axis directed 
horizontally along the central axis of the detector, such 
that beam neutrinos have a forward-going z-direction. 
The directions of the scintillator strips define a pair of 
diagonal axes U — (x + y)/y/2 and V = {—x + y)/y/2. 
Each strip provides a 2D spatial point in either the U — z 
or V — z coordinate systems, denoted the U and V views, 
respectively. 

The vertical alignment of the planes presents a source 
of difficulty in separating contained-vertex atmospheric 
neutrinos from the cosmic-ray muon background. Steep 
cosmic-ray muons incident on the detector between two 
planes can travel a significant distance into the detector 
before entering the scintillator, and therefore appear as 
contained-vertex atmospheric neutrino events. To reduce 
the background, a scintillator veto shield has been con- 
structed above the detector, and is used to tag cosmic- 
ray muons entering the detector. The veto scintillator 
modules are grouped into four sections, two per super- 
module, with a double layer on the top surface of the 
detector, and single layers diagonally above and at each 
side of the detector. To prevent gaps, adjacent modules 
overlap each other. The majority of cosmic-ray muons 
pass through two layers of scintillator before entering the 
detector and can therefore be vetoed with high efficiency. 

The veto shield is used to reject cosmic-ray muon back- 
ground in the selection of contained-vertex tracks and 
showers. An event is rejected if any activity is observed 
in the section of shield above the event vertex within a 
time window of ±50 ns. The shield efficiency is deter- 
mined using samples of cosmic-ray muons from across 
the entire data set. The fraction of cosmic-ray muons ve- 
toed by the shield is measured to be 96.6%±0.3% (syst.). 
The systematic uncertainty is obtained by modifying the 
criteria used to select the cosmic-ray muon samples and 
calculating the resulting variation in shield efficiency. 

At each stage of the contained-vertex event selection, 
the cosmic-ray muon background predictions are derived 
directly from the data by scaling down the observed dis- 
tributions of vetoed events according to the measured 
shield efficiency. A small fraction of the atmospheric 
neutrino signal is also vetoed as a result of accidental 
coincidence with noise in the shield. The loss of signal 
is determined to be 1.0%±0.2% (syst.), found by overlay- 
ing samples of veto shield data on simulated atmospheric 
neutrino events. The systematic uncertainty reflects the 
time-dependent variations in the veto shield data rates. 

Since March 2005, the MINOS Far Detector has been 
used to study neutrino interactions from the Fermilab 
NuMI accelerator beam. The beam neutrinos are iden- 
tified in the data by searching in 100 fis time windows, 



extrapolated from the beam spill times. These windows, 
which correspond to approximately 0.01% of the Far De- 
tector live time, are removed from the atmospheric neu- 
trino analysis. For the majority of running, the detec- 
tor has been magnetized to focus forward-going nega- 
tively charged muons from beam neutrino interactions. 
However, the magnetic field was reversed during the V^- 
enhanced beam running, and also for a period of several 
months prior to beam start-up, for the purpose of study- 
ing the cosmic-ray muon charge ratio [24| . 

Only data collected with both the main detector and 
veto shield fully operational are used in the analysis. The 
final data set corresponds to 2553 live-days, an exposure 
of 37.9 kton-years. 



III. ATMOSPHERIC NEUTRINO MONTE 
CARLO SIMULATION 

The MINOS Monte Carlo simulation uses separate 
programs to generate contained-vertex atmospheric neu- 
trino interactions inside the Far Detector and neutrino- 
induced muons from interactions in the surrounding rock. 
For contained-vertex atmospheric neutrino interactions, 
the NEUGEN3 simulation [25| is used to generate the in- 
teractions and hadronic final states. The transport of 
hadronic particles is then modelled using the GCALOR 
simulation 26]. For neutrino-induced rock- muons, the 
nuance generator [27|_is used, with the GRV94 parton 
distribution functions 28| . The NUANCE simulation then 
propagates the muons from the rock to the edges of the 
detector. For both Monte Carlo samples, a GEANT3 [29[ 
simulation of the Far Detector is used to model particle 
transport and detector response. 

For contained-vertex atmospheric neutrino events, the 
simulation uses the flux calculation of Barr et al. (30j 
(Bartol 3D). For neutrino energies below 10 GeV, this cal- 
culation is based on a 3D simulation, with separate flux 
tables provided for the cases of solar minimum and solar 
maximum. Above 10 GeV, a ID simulation has been used 
and a single set of flux tables is provided. For neutrino- 
induced rock-muons, where the majority of parent neu- 
trinos have energies greater than 10 GeV, the simulation 
uses the earlier ID calculation by the Bartol group [3l[ 
(Bartol ID). The rock-muons are then reweighted as a 
function of their parent neutrino energy and zenith angle 
using the ratio of the Bartol 3D and Bartol ID fluxes. 
Both the contained-vertex and rock-muon samples use 
distributions of neutrino production height obtained from 
a separate simulation of cosmic-ray interactions in the at- 
mosphere and parameterized in terms of neutrino energy 
and zenith angle [5]. 

The Far Detector data set spans a significant fraction 
of the full solar cycle. For atmospheric neutrinos above 
500 MeV (the MINOS energy threshold), the predicted 
neutrino interaction rate is 7% higher at solar minimum 
than solar maximum. These solar cycle effects are ac- 
counted for by taking a weighted average of the Bartol 3D 
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FIG. 1. The simulated atmospheric neutrino energy spec- 
trum in the MINOS Far Detector for 2553 live-days of data. 
Separate distributions are plotted for contained-vertex neu- 
trino interactions and neutrino-induced rock-muons showing 
the predictions for the case of no oscillations, and for os- 
cillations with Am 2 = 2.32 x 10" 3 eV 2 and sin 2 20 = 1.0. 
The contained-vertex neutrino interactions are generated in 
the range 0.2 — 50 GeV, with a median value of 2 GeV; the 
neutrino-induced rock-muons range up to neutrino energies of 
10 TeV, with a median value of 50 GeV. The effect of — !> v T 
oscillations is visible for neutrino energies below 100 GeV. 



fluxes calculated at solar minimum and solar maximum. 
The variation of solar activity over time is determined 
by parameterizing available atmospheric neutron data 
from the CLIMAX experiment [321. Combining these 
data with the Far Detector running periods, the fluxes 
at solar minimum and solar maximum are combined in 
proportions 70% and 30%, respectively. 

Figure [T] shows the simulated atmospheric neutrino 
energy spectrum for contained-vertex neutrino interac- 
tions and neutrino-induced rock-muons, plotted for the 
case of no oscillations, and for oscillation parameters of 
Am 2 = 2.32 x 10" 3 eV 2 and sin 2 20 = 1.0 12]. These os- 
cillation parameters are also used to calculate predicted 
atmospheric neutrino event rates in Sections IV1 and PVtl of 
this paper. 



A. Systematic uncertainties in atmospheric 
neutrino simulation 

The predicted atmospheric neutrino event rates have 
large uncertainties arising from the atmospheric neutrino 
flux and interaction models. The Bartol group has car- 
ried out a detailed study of the systematic uncertainties 
in their 3D flux model |33| . These flux uncertainties have 
also been studied by comparing the Bartol model with 
the alternative 3D calculations of Battistoni et al. (34| 
and Honda et al. (35j . The systematic uncertainties used 



in this analysis are based on the results of the Bartol 
study but are also found to cover the differences between 
the different flux models. 

The dominant source of uncertainty in the overall rate 
of contained-vertex neutrinos and neutrino- induced rock- 
muons is the systematic uncertainty in the normaliza- 
tion of the atmospheric neutrino flux simulation. This 
overall uncertainty increases with neutrino energy due to 
rising uncertainties in the primary cosmic-ray flux and 
hadroproduction models. For contained-vertex neutri- 
nos, which have a median energy of 2 GeV and lie primar- 
ily below 10 GeV, an overall uncertainty of 15% is applied 
in this analysis. For neutrino-induced rock-muons, where 
the parent neutrino has a median energy of 50 GeV, with 
an energy spectrum that ranges up to 10 TeV, a larger 
uncertainty of 25% is applied. 

At low neutrino energies, many of the systematic un- 
certainties in the flux model cancel in the ratios of dif- 
ferent flux components. In the 1 — 5 GeV region, the 
uncertainties in the {vp,+v^)/(y e -\-T> e ) and v^jv^ flux ra- 
tios, and in the up-down ratio of upward-going neutri- 
nos to downward-going neutrinos, are calculated to be 
smaller than 5% [3J|. For the analysis presented here, 
the uncertainty in the v^jv^ ratio is of greatest impor- 
tance. At energies below 10 GeV, cosmic-ray hadropro- 
duction predominantly yields pions, with each charged 
pion producing a single pair of and in its decay 
chain. Therefore, the v^jv ^ ratio approaches unity, with 
a high degree of cancellation in its systematic uncertainty. 
For the contained-vertex atmospheric neutrino sample, 
a conservative uncertainty of 4% is placed on this ra- 
tio. At energies above 10 GeV, the cancellations in the 
ratio quickly diminish, as the large uncertainty in the 
kaon component of hadroproduction becomes a signifi- 
cant factor, and an increasing fraction of atmospheric 
muons strike the ground before decaying. Therefore, for 
the neutrino-induced rock-muon sample, a larger uncer- 
tainty of 10% is placed on the v^/v^ ratio. 

Additional systematic uncertainties in the predicted 
atmospheric neutrino event rate arise from the neutrino 
interaction model. The uncertainty in the total CC 
cross-section peaks at 8% in the 1 — 5 GeV region [3(| , 
corresponding to the transition region between the mod- 
els of quasi-elastic and resonance neutrino interactions. 
At higher energies, where deep inelastic interactions are 
dominant, the total interaction cross-section is well con- 
strained by experiment, and the uncertainty falls to 2%. 
For antineutrinos, where there is limited experimental 
data below 5 GeV, the predicted CC cross-section has 
a larger uncertainty. For the analysis presented here, 
an energy-dependent systematic uncertainty band on the 
Vp-I^n cross-section ratio has been calculated by varying 
the input parameters to the NEUGEN3 interaction model 
according to their given uncertainties [3?} ■ The aver- 
age uncertainty in the v^jv ^ cross-section ratio is then 
calculated by integrating across this uncertainty band, 
weighting each bin of neutrino energy by the predicted 
rate of atmospheric neutrino and v ^ CC interactions. 
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For the contained-vertex neutrino sample, this procedure 
yields an overall uncertainty of 8.5%. For the neutrino- 
induced rock-muon sample, where the majority of events 
are produced by deep inelastic neutrino interactions, the 
calculation returns a smaller uncertainty of 4%. 

The atmospheric neutrino event rate at the Soudan 
mine has been previously measured by the Soudan 2 ex- 
periment. The Soudan 2 analysis of atmospheric elec- 
tron neutrinos indicates that the predicted interaction 
rates obtained by combining the Bartol 3D flux model 
and NEUGEN cross-section model should be scaled by 
0.91 ±0.07 @. However, although the Soudan 2 and MI- 
NOS detectors are located at the same site, Soudan 2 has 
a lower neutrino energy threshold of 300 MeV, compared 
with 500 MeV for MINOS. An analysis of contained- 
vertex showers from atmospheric neutrinos by MINOS, 
based on 418 live-days of data, yields a scale factor of 
1.08 ± 0.12 (stat.) ± 0.08 (syst.) [3j|. The systematic un- 
certainties applied in this analysis cover both these mea- 
surements. 



IV. EVENT RECONSTRUCTION 

The data are reconstructed using an algorithm which 
identifies the track and shower topologies in each 
event [39]. Reconstructed tracks typically contain hits 
in one strip per plane and are principally produced by 
muons; reconstructed showers contain hits in multiple 
strips per plane and are produced by hadronic and elec- 
tromagnetic particles. 

Initially, the particle tracks and showers are recon- 
structed independently in each of the U and V views; 
these 2D views are then matched to generate a 3D event. 
A Kalman filter algorithm is used to determine the tra- 
jectory of each muon track, accounting for energy loss 
in the detector and curvature in the magnetic field [40l |. 
This algorithm also reconstructs the start and end points 
of each track, which are distinguished using timing infor- 
mation. For tracks where the end point lies inside the 
fiducial volume, the muon momentum is reconstructed 
from the measured track length; for exiting tracks, the 
momentum is obtained from the fitted track curvature. 
In both cases, the fitted curvature is used to determine 
the muon charge sign. For atmospheric neutrino events 
containing a reconstructed track, the interaction vertex 
is given by the start point of the track; if there is only a 
reconstructed shower, the vertex is given by the centroid 
of the shower. 

The propagation direction of each muon along its re- 
constructed track is determined using timing informa- 
tion. The MINOS Far Detector has a single-hit tim- 
ing resolution of approximately 2.5 ns, which enables the 
muon direction to be reconstructed with high purity for 
tracks spanning ten or more scintillator planes. The Far 
Detector timing system is calibrated using cosmic-ray 
muons, which are used to determine the time offsets in 
each readout channel and to correct for shifts in these 



offsets resulting from swapped readout components (39|. 

The detector is calibrated using a combination of LED 
light injection and the average pulse height response of 
each strip using cosmic ray muons [23]. A minimum- 
ionizing muon passing through a scintillator strip at 
normal incidence generates a combined signal of ap- 
proximately 10 photo-electrons (PEs). The selection of 
contained-vertex tracks and showers makes use of the en- 
ergy profile of events. 

For contained-vertex atmospheric and events, 
the emitted muon is typically accompanied by some re- 
constructed shower activity at the interaction vertex, 
produced by the hadronic system. The total hadronic 
energy is determined by summing the calibrated pulse 
heights in the reconstructed shower. For low energy 
showers, large fluctuations can occur, degrading the 
hadronic energy resolution. To reduce the size of these 
fluctuations, the pulse heights are first raised to a power 
before being summed together. The exponent used in 
this procedure is increased as a function of shower en- 
ergy from a minimum of 0.25 at the lowest shower en- 
ergies to a maximum of 1.0 for shower energies above 
18GeV [4l[. Studies of simulated atmospheric neutrinos 
show that, relative to a linear summation of pulse heights, 
the hadronic energy resolution improves from 55% to 45% 
for reconstructed showers in the 1 GeV region. 



V. EVENT SELECTION 

An initial selection is applied to all events, ensuring 
a good reconstruction quality. The selected events are 
then separated into a track-like sample containing re- 
constructed tracks that span 8 or more planes, and a 
shower-like sample containing reconstructed showers that 
span 4 or more planes. The track-like sample is used 
for the selection of contained-vertex muons and neutrino- 
induced rock-muons; the shower-like sample is used for 
the selection of contained-vertex showers. Initially, 2% 
of events are placed in both the track-like and shower- 
like samples. Any duplicate events are removed from the 
shower-like sample after the full selection has been ap- 
plied. After the initial selection, the observed event rate 
is 55,000 events/day, dominated by cosmic-ray muons. 
The corresponding predicted atmospheric neutrino rates 
are 0.8 events/day from contained-vertex interactions, 
and 0.3 events/day from neutrino-induced rock-muons, 
after accounting for oscillations. 

The atmospheric neutrino signal is separated from the 
cosmic-ray background using two characteristic signa- 
tures of atmospheric neutrino interactions: either a re- 
constructed vertex inside the fiducial volume; or a recon- 
structed upward-going or horizontal muon trajectory. A 
set of requirements on event containment and topology is 
first applied to select contained-vertex tracks and show- 
ers, using the veto shield to reduce the cosmic-ray muon 
background. A set of requirements on event timing infor- 
mation and length is then applied to select upward-going 
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FIG. 2. Distributions of the trace variable, Az, for contained- 
vertex tracks. This estimates the distance in z travelled by 
a cosmic-ray muon inside the detector before first entering 
the scintillator. The hatched histogram shows the simulated 
prediction for the atmospheric neutrino signal; the solid line 
shows the predicted total rate, given by the sum of the sig- 
nal and the cosmic-ray muon background; the points show 
the observed data. The background distribution is peaked to- 
wards low values of Az, and the arrow indicates the selection 
applied to reduce the background. 

and horizontal muons produced by neutrino interactions 
in the detector or surrounding rock. The full selection is 
described in the following sections. 

A. Selection of contained-vertex tracks 

The contained-vertex track selection criteria identify 
and F M CC atmospheric neutrinos. For this sample, 
the veto shield selection is first used to reduce the level 
of cosmic-ray muon background. A set of containment 
and topology selection criteria are then applied to the 
remaining tracks. For each track, the reconstructed tra- 
jectory has two ends, corresponding to the first and last 
scintillator hits on the track. Since cosmic-ray muons 
are incident from above, the majority of these selection 
criteria are applied at the upper end of the track. The 
following selection criteria are applied [39|, [42|, EH : 

1. Fiducial Cuts: The reconstructed track vertex is re- 
quired to lie within a fiducial volume starting 0.2 m 
inside any edge of the detector, 5 planes from the 
ends of each super-module, and 0.4 m from the cen- 
ter of the coil hole. To reject cosmic-ray muons 
that enter the detector through the coil hole, the 
coil cut is increased to 1 m in the first and last 20 
planes of the detector. In addition, if the vertex 
is reconstructed at the lower end of the track from 
timing information, these selection criteria are also 
applied at the upper end of the track. 



2. Trace Cut: The cosmic-ray muons that pass the 
fiducial requirements typically enter the detector 
at a small angle to the planes, and can travel a 
significant distance through the detector before en- 
tering the scintillator. However, the distance trav- 
elled along the z-axis is typically small, and so the 
background can be reduced by placing a minimum 
requirement on this distance. For each event, a 
detector entry point is estimated by extending the 
reconstructed trajectory at the upper end of the 
track upwards to the edge of the detector. The dis- 
placement along the z-axis, Az, between the entry 
point and upper end of the track is then calculated 
(this quantity is referred to as the "trace"). Fig- 
ure [5] shows the predicted and observed distribu- 
tions of the trace variable. The cosmic-ray muon 
background is peaked towards low values, whereas 
the atmospheric neutrino signal distribution has 
a flatter distribution. To reduce the background, 
events are required to satisfy Az > 0.5 m. 

3. Topology Cuts: The cosmic-ray muon background 
events that pass the trace cut typically travel a 
significant distance in a single steel plane and its 
associated air gap before entering the scintillator. 
A number of these cosmic-ray muons undergo sig- 
nificant bending in the magnetic field, and some 
muons reverse the horizontal component of their 
direction. As a result, the reconstruction may miss 
the first hit on the track or underestimate the steep- 
ness of the track. These events are characterized 
by clusters of hits above the upper end of the re- 
constructed track. To reduce this background, the 
charge-weighted mean and RMS displacements of 
strips, denoted (Ajjv) and (Afjy)?, are calculated 
separately in the U and V views for a ±4 plane 
window around the upper end of the track. Events 
are rejected if (Ajjv) > 0.25 m, indicating that sig- 
nificant energy has been deposited above the track, 
or if (Afjy)? > 0.5 m, indicating that there was 
significant scatter at the upper end of the track. 
A set of 3D displacements is also calculated using 
the same ±4 plane window, by combining all possi- 
ble pairs of U and V strips in adjacent scintillator 
planes. The maximum 3D displacement from the 
upper end of the track, A 1 ^ ax , is then calculated, 
and events are rejected if AJj: ax > 1.25 m. 

4. Pulse Height and Direction Cuts: The cosmic-ray 
muon background is also characterized by large de- 
posits of energy at the upper end of the track, due 
to the long distance travelled in the first plane. 
This background is reduced by finding the max- 
imum pulse height, Q v t x , in a ±4 plane window 
around the upper end of the track. Events are re- 
jected if Q v tx > 300 PEs. The Q v t x requirement 
is tightened to 75 PEs if the track is both short, 
spanning fewer than 25 planes, and steep, satisfying 
cos 8 y > 0.7 or | cos 02 1 < 0.5. Here, 9 y and 8 Z are 
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FIG. 3. Distributions of the pulse height at the upper end of the track (Q v tx), plotted against the z-component (cos# z ) and 
y-component (cos8 y ) of the downward track direction. The distributions are plotted for contained- vertex muons that pass the 
trace and topology requirements and span fewer than 25 planes. The plots on the left show simulated atmospheric neutrinos; 
those on the right show the cosmic-ray muon background. The background events are associated with large pulse heights and 
directions parallel to the vertically-aligned scintillator planes. The hatched area is the region rejected by the pulse height and 
direction selection criteria as part of the topology requirements. Note that the requirement of a track spanning > 8 planes 
causes the acceptance to drop to zero as | cos# z | approaches 0, and as cos8 y approaches —1. 



taken as the angles between the reconstructed tra- 
jectory at the upper end of the track, and the y and 
z axes, respectively. Figure [3] shows the predicted 
signal and background distributions of Q v t x as a 
function of 8 y and 9 Z for the short tracks, indicat- 
ing the selection requirements applied to separate 
the atmospheric neutrino signal from the cosmic- 
ray muon background. 

The track containment requirements yield 801 events 
from the data. This compares with total predictions 
of 934 ± 134 events for no oscillations, and 698 ± 99 
events for oscillations with Am 2 = 2.32 x 10~ 3 eV 2 and 
sin 2 2# = 1.0. The uncertainties in these predictions are 
dominated by the 15% uncertainty in the overall normal- 
ization of the contained-vertex Monte Carlo simulation. 
The combined and T>^ CC components form 92% of 
the total predicted event rate before oscillations. This 
component oscillates in the two-flavor model and there- 
fore represents the signal in the oscillation analysis. The 
combined i> e +V e CC and NC components, which do not 
oscillate in the two-flavor model, form a 5% background. 
The cosmic-ray muon prediction of 34 ± 3 events corre- 



sponds to a 3% background level. 

B. Selection of upward-going and horizontal tracks 

For upward and horizontal angles, where the rock over- 
burden exceeds 14,000 m water-equivalent, the absorp- 
tion of cosmic-ray muons by the earth is sufficiently high 
that the observed flux of muons is dominated by atmo- 
spheric muon neutrino interactions [44 |. At the Soudan 
mine, this corresponds to zenith angles in the range 
cos6 z < 0.14 [45]]. Therefore, upward-going and horizon- 
tal muons in the MINOS Far Detector provide a signature 
for atmospheric neutrinos. 

Upward-going and horizontal muons are selected based 
on the reconstructed zenith angle at the track vertex. 
The direction of muon propagation along the track is re- 
constructed using timing information. This is then used 
to distinguish between the track vertex and end points. 
The reconstructed track vertex can either be inside or 
outside the fiducial volume. Hence, this sample of events 
provides a source of both contained-vertex muons and 
neutrino-induced rock-muons. 
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FIG. 4. Distribution of 1//3 normalized velocity variable, 
demonstrating the purity of the track direction identification. 
The 1//3 variable is the gradient of a linear fit to the measured 
times as a function of distance along each track. The distribu- 
tion is plotted for all tracks that pass the topology and timing 
selections. The peak at —1.0 corresponds to downward-going 
muons; the peak at +1.0 to upward-going muons. A good 
separation is achieved between the upward-going neutrino- 
induced signal, and downward-going cosmic-ray muon back- 
ground. 

To determine the track direction from timing informa- 
tion, two linear fits are applied to the measured times of 
the track hits, as a function of their distance along the 
track. The gradients are constrained to be ±l/c, corre- 
sponding to forward and backward propagation along the 
track at the speed of light, c. The hits are weighted as 
a function of pulse height to account for the variation in 
the single-hit timing resolution, which is better for larger 
pulse heights due to increased photon statistics. For each 
propagation direction, the goodness of the timing fit is 
given by its RMS timing residual. The smaller of the two 
RMS values is labelled r^, and the larger is labelled Th- 
The propagation direction is determined by the timing 
fit with the smaller RMS residual. 

The neutrino-induced muons must be separated from 
a high background of cosmic-ray muons whose direction 
is mis-reconstructed. To ensure that the track direction 
is reconstructed unambiguously, the following selection 
criteria are applied: 

1. Topology Cuts: To ensure that there are suffi- 
cient hits on the track for the zenith angle to be 
reconstructed accurately, and the propagation di- 
rection determined unambiguously, reconstructed 
tracks are required to span more than 15 planes 
and to travel more than 1.5 m. For upward- going 
tracks, the track vertex point is required to lie be- 
low the track end point, or, to account for possible 
track curvature, no more than 0.5 m above it. For 




cos e z 

FIG. 5. Distribution of reconstructed zenith angle for muons 
with good timing and topology. In the range cos# z > 0.10, 
the observed rate of muons is dominated by the cosmic-ray 
background and falls steeply as the mean rock overburden in- 
creases rapidly. For cos6 l z < 0.10, the distribution flattens, as 
the cosmic-ray muon flux falls below that of neutrino-induced 
muons. To minimize the background from cosmic-ray muons, 
events are required to satisfy: cos# z < 0.05. 

downward-going tracks, the track vertex require- 
ment is reversed. 

2. Timing Cuts: To ensure that the muon propa- 
gation direction is identified unambiguously from 
timing information, a set of selection requirements 
are placed on the quality of the linear timing fits. 
The difference between the two RMS residuals is re- 
quired to satisfy tl ~ r H < —1.66 ns, a significant 
fraction of the detector timing resolution. In addi- 
tion, an upper requirement of tl < 4.66 ns is placed 
on the smaller residual and a lower requirement of 
Th > 3.66 ns is placed on the larger residual. The 
ratio between the best fit RMS residual, n,, and 
the track length, I, also provides a means of select- 
ing events with well-measured timing information. 
Events are required to satisfy r^filjc) < 0.577. 

These selection criteria identify clean samples of both 
upward-going and downward-going muons from the data. 
As a check on the quality of separation between these 
two samples, an additional unconstrained linear timing 
fit is applied to the selected events. The measured times 
along the track are fitted as a function of their upward 
distance along the track. The fits return a gradient, 1/v, 
where v is the reconstructed velocity. Figure 2] shows the 
distribution of the normalized gradient, 1/(3 = l/(v/c). 
A good separation is achieved between the upward-going 
events, dominated by neutrino-induced muons, and the 
downward-going events, dominated by the cosmic-ray 
muon background. 

Figure [5] shows the distribution of reconstructed zenith 
angle for the selected events. In the region cos 6^ > 0.10, 
the event rate falls steeply with zenith angle, as the 
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rapidly increasing rock overburden reduces the incident 
flux of cosmic-ray muons. In the region cos9 z < 0.10, 
the event rate flattens and becomes approximately con- 
stant, as neutrino-induced muons become the dominant 
flux component. A residual background arises from low 
momentum cosmic-ray muons that deflect significantly 
due to multiple Coulomb scattering in the rock and en- 
ter the detector in a horizontal direction. To minimize 
the cosmic-ray muon background, the selected sample of 
upward-going and horizontal muons is required to satisfy: 
cos0 z < 0.05. An estimate of the remaining background 
is obtained from an exponential fit to observed data in 
the region | cos(9 z | < 0.20. This method returns a back- 
ground prediction of 0.5 selected events. The cosmic-ray 
muon background component is neglected in the subse- 
quent analysis. 

Overall, 665 upward-going and horizontal muons are 
selected from the data. This compares with total predic- 
tions of 882 ± 146 events in the absence of oscillations, 
and 623 ±113 events for input oscillation parameters of 
Am 2 = 2.32 x 10" 3 eV 2 and sin 2 26> = 1.0, where the 
uncertainties are dominated by the normalizations of the 
contained- vertex and neutrino-induced rock-muon Monte 
Carlo simulations. The selected sample divides into 466 
neutrino-induced rock-muons and 199 contained-vertex 
muons. The latter sample contains 95 events already 
selected as contained-vertex muons by the containment 
requirements described above, along with an additional 
104 upward-going and horizontal events. 

C. Selection of contained-vertex showers 

The contained-vertex shower sample is primarily com- 
posed of the v e +V e CC and NC atmospheric neutrino 
component. For this sample, the main background arises 
from cosmic-ray muons incident at steep angles, which 
radiate large showers and span a small number of planes. 
The veto shield is first applied to reduce cosmic-ray muon 
background. The following selection criteria are then ap- 
plied to the remaining showers (38[: 

1. Fiducial and Trace Cuts: The fiducial require- 
ments described above for contained-vertex tracks 
are also applied in the selection of contained-vertex 
showers. For reconstructed showers, the vertex res- 
olution is poorer and cosmic-ray background level 
larger. Therefore, a tighter fiducial requirement of 
0.4 m is applied at each edge of the detector for this 
sample. The shower direction also has a poorer res- 
olution, and therefore a tighter trace cut of 0.8 m is 
applied. Figure [6] shows the predicted and observed 
Az distributions for contained-vertex showers. 

2. Topology Cuts: Further selection requirements are 
applied to identify the characteristic topologies of 
neutrino-induced showers. The longitudinal profile 
of these showers typically rises to a maximum and 
then falls smoothly, whereas showers generated by 
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FIG. 6. Distributions of the trace variable, Az, for contained- 
vertex showers. This variable estimates the distance in z trav- 
elled by a cosmic-ray muon inside the detector before entering 
the scintillator. The hatched histogram shows the simulated 
prediction for the atmospheric neutrino signal. The solid line 
gives the predicted total rate, dominated by the cosmic-ray 
muon background. The points show the observed data. The 
background is peaked towards low values of Az , since cosmic- 
ray muons typically travel a small distance in z before entering 
the scintillator. The arrow indicates the selection applied on 
Az to reduce the background. 

cosmic-ray muons typically deposit a lot of energy 
in a single plane, or contain large fluctuations be- 
tween planes. The cosmic-ray background is also 
found to be higher for shorter showers. Therefore, 
events are separated into short (< 8 planes) and 
long (> 8 planes) samples, and tighter selection 
criteria are applied to the short sample. To char- 
acterize the shower topology, the mean and RMS 
number of strips per plane {{Wuv) and (W^ v )^), 
and pulse height per plane ((Q s hw) and (Q 2 shw )^), 
are calculated. These shower topology variables are 
required to satisfy: 

(Wuv) < 5 (4) strips, 
(W* v )? < 4 (3) strips, 
{Q.hw) < 150 (100) PEs, 
(Q 2 ,J3 < 150 (100) PEs, 

for long (short) showers, respectively. Figure [7] 
shows the predicted and observed distributions of 
the mean and RMS shower pulse height variables. 
An analysis of the principal moments of the shower 
is also used to distinguish the shower-like event 
topology of atmospheric neutrinos from the track- 
like topology of cosmic-ray muons. The moment 
of inertia tensor is constructed from the relative 
positions of the shower strips, weighted by their 
pulse height. The tensor is diagonalized and selec- 
tion criteria are placed on the largest eigenvalue, 
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FIG. 7. Distributions of mean and RMS shower pulse height per scintillator plane observed for contained-vertex showers that 
pass the fiducial and trace cuts. The hatched histogram shows the simulated prediction for the atmospheric neutrino signal. 
The solid line gives the total prediction, dominated by the cosmic-ray muon background. The points show the observed data. 
The arrows indicate the selection applied to reduce the background. 
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FIG. 8. Distributions of reconstructed neutrino energy, plot- 
ted for selected contained-vertex showers. The dotted line 
shows the prediction for no oscillations; the solid line shows 
the prediction for Am 2 = 2.32 x 10~ 3 eV 2 and sin 2 26» = 1.0; 
the shaded histogram shows the predicted cosmic-ray muon 
background; the points with errors show the observed data. 
The v^+V^ CC component is small and hence the total predic- 
tion does not depend strongly on the oscillation parameters. 



3. Removal of Selected Tracks: After applying the 
above selection criteria, it is found that 2% of the 
resulting events have been previously selected as 
contained-vertex tracks. Monte Carlo studies indi- 
cate that approximately half of the duplicate events 
are or CC interactions. As a hnal step, 
these events are removed from the contained-vertex 
shower sample. 



The shower containment requirements select 701 
events from the data, compared with a total predic- 
tion of 727 ± 101 events in the absence of oscillations, 
and a prediction of 684 ± 95 events for oscillations with 
Am 2 = 2.32 x 10~ 3 cV 2 and sin 2 26< = 1.0. The uncertain- 
ties are dominated by the 15% uncertainty in the overall 
normalization, but also includes additional uncertainties 
of 20% and 5% in the NC and v e +v e CC components, 
respectively. The cosmic-ray muon prediction of 87 ± 9 
events represents a background level of 12%. Figure [5] 
shows the predicted and observed energy distributions of 
the selected events. Since the selected sample contains a 
small v^+Vu CC component, the distribution has only a 
weak dependence on the oscillation parameters. 
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D. Summary of results from atmospheric neutrino 
event selection 

In total, 2072 candidate atmospheric neutrino events 
are selected from the data. For analysis, the events are 
grouped into: 905 contained-vertex muons, with vertex 
positions inside the fiducial volume (including both fully- 
contained and partially-contained muons); 466 neutrino- 
induced rock-muons, with vertex positions outside the 
fiducial volume; and the 701 contained-vertex showers. 
Table Q] gives the predicted event rates for each of these 
samples. 



E. Selection of high resolution event sample 

A sample of high resolution contained-vertex muons, 
with well-measured muon propagation direction, is se- 
lected from the data. These high resolution events are 
required to satisfy minimum track length requirements 
of 10 planes and lm. They are also required to pass a 
timing requirement of — rjy < —0.66 ns. These selec- 
tion criteria are found to correctly distinguish the track 
direction in 99% of simulated atmospheric neutrinos. 

A total of 631 high resolution contained-vertex muons 
are selected, with the remaining 274 contained-vertex 
muons classified as low resolution. In the high resolu- 
tion sample, 261 events are classified as upward-going 
and 370 events are classified as downward-going. The 
measured up-down ratio is i?^ a = 0.71 ± 0.06 (stat.), 
where the statistical error corresponds to the 68% con- 
fidence interval calculated using Poisson statistics [46| . 
The predicted ratio, calculated from the simulation, in 
the absence of oscillations, is R^S = 1-14 ± 0.03 (syst.). 
The 3% systematic uncertainty combines the uncertain- 
ties in the event selection, and the atmospheric neutrino 
flux simulation. The double ratio between the observed 
and predicted up-down ratio is: 

Ru/d/Ru/d = 0-62 ± 0.05 (stat.) ± 0.02 (syst.). 

This ratio is in excess of 6 standard deviations from unity, 
indicating the presence of neutrino oscillations. 



VI. SEPARATION OF NEUTRINOS AND 
ANTINEUTRINOS 

The high resolution contained-vertex muon sample and 
neutrino-induced rock-muon sample are separated into 
candidate neutrinos and antineutrinos based on the re- 
constructed muon charge sign. The Kalman filter returns 
a best fit value of q/p, and its uncertainty, <J q / p , where q 
is the muon charge sign and p is the muon momentum. 
The selected events are classified as neutrinos if q < 0, 
and as antineutrinos if q > 0. 



Two criteria are used to select events with significant 
track curvature and therefore well-measured charge sign. 
First, a requirement is placed on the relative size of the 
track fit uncertainty, \q/p\/ <J q / pi which indicates the sig- 
nificance of the track curvature @ . Events are required 
to satisfy \q/p\/cr q / p > 2.5. The reconstructed track is 
then used to calculate a variable measuring the straight- 
ness of the track, with the aim of excluding tracks which 
do not have significant curvature. A straight line is drawn 
between the reconstructed start and end point of the 
track, and a chi-squared variable, xfme/^ /' ^ s calculated 
from the deviations of the track strips from this line (tJ. 
Events are required to satisfy xjine/dof > 4.0. Figure [9] 
shows the observed and predicted distributions of these 
two selection variables for contained-vertex muons and 
neutrino-induced rock-muons. 

The charge selection criteria are found to correctly 
identify the muon charge in 97% of simulated contained- 
vertex interactions, and 99% of simulated neutrino- 
induced rock-muons. The selection efficiencies are 87% 
and 59%, respectively, calculated as a fraction of the 
number of selected events with well-measured direction. 
The lower efficiency for neutrino-induced rock-muons re- 
flects their higher average momentum, resulting in more 
events with ambiguous track curvature. 

Table |TT] gives the predicted and observed numbers of 
neutrinos and antineutrinos for each category of event. 
For contained-vertex muons, the charge-separation pro- 
cedure returns 379 neutrinos and 173 antineutrinos, giv- 
ing a measured charge ratio of R^ 1 ® = 0.46^°'q4 (stat.). 
The predicted value of the charge ratio is calculated from 
the simulation to be R¥P = 0.49 ±0.05 (syst.). The pre- 
diction is almost entirely independent of any input os- 
cillations, provided that equal parameters are used for 
neutrinos and antineutrinos. The overall systematic un- 
certainty of 10% is obtained by combining uncertainties 
of 8.5% in the ratio of the neutrino and antineutrino 
interaction cross-sections, 4% in the flux ratio of neu- 
trinos and antineutrinos, and 3% in the purity of the 
charge separation. The double ratio between the ob- 
served and predicted charge ratios is calculated to be: 
R±f*/R¥C = 0.93 ± 0.09 (stat.) ± 0.09 (syst.). 

For neutrino-induced rock-muons, 152 neutrinos and 
95 antineutrinos are selected, giving a measured charge 
ratio of Rif^ = 0.63^;^ (stat.). The predicted value of 
the charge ratio is calculated from the simulation to be 
= 0.48 ±0.06 (syst.). The overall systematic uncer- 
tainty of 12.5% is obtained by combining uncertainties 
of 10% in the flux ratio, 4% in the cross-section ratio, 
and 6% in the charge-separation purity. The double ra- 
tio between the observed and predicted charge ratio is: 
Ri% a /R^/ C v = l-29±8:i? (stat.) ± 0.16 (syst.). 

The charge-separated samples of contained-vertex 
muons and neutrino-induced rock-muons are combined 
to give an overall double ratio of: 

R±f v a /R^c = 1.03 ± 0.08 (stat.) ± 0.08 (syst.). 
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Data 






Prediction (no oscillations) 










Cosmic-ray /j, 


^V+^m cc 


v e +v e CC v T +v T CC 


NC 


Rock-/x 


Total 


Contained-vertex muons 


905 


34 ± 3 


998 ± 150 


35 ±6 


25 ±6 


9 ± 2 


1100 ± 159 


Neutrino-induced rock-muons 


400 




26 ±4 


0±0 


0±0 


£ A A -J_ 1QC 

044 ± loD 


f U ± loo 


Contained-vertex showers 


701 


87 ± 9 


157 ± 24 


358 ± 57 


124 ±31 


1 ± 


727 ± 101 


Total 


2072 






2397 ± 296 










Data 




Prediction (Am 2 = 2.32 x 10" a eV 


2 , sin 2 26> = 


:1.0) 








Cosmic-ray 


^ M + !7 M CC 


v e +v e CC v T +v T CC 


NC 


Rock-/i 


Total 


Contained-vertex muons 


905 


34 ±3 


689 ± 103 


35 ± 6 3 ± 1 


25 ±6 


6± 1 


792 ± 113 


Neutrino-induced rock-muons 


466 




14 ±2 


0±0 0±0 


0±0 


433 ± 108 


447 ± 108 


Contained-vertex showers 


701 


87 ±9 


110 ± 16 


358 ± 57 5 ± 2 


124 ±31 


0±0 


684 ± 95 



Total 2072 1923 ± 235 



TABLE I. Summary of atmospheric neutrino selection, separated into the different categories of selected and simulated event. 
The Monte Carlo predictions are given separately for contained-vertex atmospheric neutrinos and neutrino-induced rock-muons, 
with the contained-vertex predictions also separated by neutrino interaction type. The cosmic muon background prediction 
is calculated directly from the data by weighting vetoed events according to the measured shield efficiency. In the top table, 
the predictions are calculated in the absence of neutrino oscillations; the bottom table uses representative oscillation parameters 
of Am 2 = 2.32 x 10" 3 eV 2 and sin 26 = 1.0. The predictions and their uncertainties have been rounded to the nearest event. 
Note that many of the uncertainties are correlated and cancel in the ratios and fits described in this paper. 





Data 






Prediction (no 


oscillations) 










i> CC 


v„ cc 


Rock-/i~ 


Rock-/i + 


Other 


Total 


Contained-vertex muons (^t~) 


379 


425 ± 64 


4± 1 


4± 1 


0±0 


13 ± 1 


445 ± 65 


Contained-vertex muons (/i + ) 


173 


12 ±2 


190 ± 28 


0±0 


1±0 


15 ±2 


219 ± 31 


Neutrino- induced rock-muons 


152 


16 ±2 


0±0 


215 ± 54 


1±0 


0±0 


233 ± 54 


Neutrino-induced rock-muons (/i + ) 


95 


0±0 


8± 1 


3± 1 


102 ± 25 


0±0 


112 ± 26 




Data 




Prediction (Am 2 = 2.32 x 10 -3 eV 2 , sin 2 


29 = 1.0) 








i/p CC 


V„ CC 


Rock-/i 


Rock-/^ + 


Other 


Total 


Contained-vertex muons (/i~) 


379 


294 ± 44 


3±0 


2± 1 


0±0 


15 ±2 


314 ± 46 


Contained-vertex muons (fJ> + ) 


173 


9± 1 


132 ± 20 


0±0 


1±0 


16 ±2 


158 ± 22 


Neutrino- induced rock-muons (/z ) 


152 


9± 1 


0±0 


151 ± 38 


1±0 


0±0 


161 ± 38 


Neutrino-induced rock-muons (fi + ) 


95 


0±0 


4± 1 


2±0 


68 ± 17 


0±0 


74 ± 18 



TABLE II. Results from the separation of contained-vertex and neutrino-induced rock-muons by reconstructed charge sign, 
into selected samples of neutrinos and antineutrinos. The predictions from each simulated sample are given separately for true 
neutrinos and antineutrinos; the column labelled 'Other' is the sum of the cosmic-ray muon, i/ e +V e , v T +V T and NC backgrounds. 
All the predictions and their uncertainties have been rounded to the nearest event. Note that many of the uncertainties are 
correlated and cancel in the ratios and fits described in this paper. 



This result is consistent with unity. 



VII. OSCILLATION ANALYSIS 

Two oscillation fits are applied to the selected data. 
The first is a two-parameter fit, which outputs equal os- 
cillation parameters for neutrinos and antineutrinos; the 
second is a four-parameter fit, which outputs separate 
oscillation parameters for neutrinos and antineutrinos. 
Each fit is applied to the reconstructed L/E distributions 
of selected neutrinos and antineutrinos. The neutrino 
propagation length, L, is determined from the recon- 
structed zenith angle of the muon track. For contained- 
vertex muons, the parent neutrino energy, E, is found 
by summing the reconstructed muon energy and visible 
shower energy. For neutrino-induced rock-muons, only 



the muon energy is used, due to the increased uncertain- 
ties in the shower simulation and calibration at the edge 
of the detector, or because the interaction occurs outside 
the detector and hence the vertex is not visible. Fig- 
ure [TU] shows the predicted and observed reconstructed 
zenith angle distributions for contained-vertex muons 
and neutrino-induced rock-muons. Figures [TT] and 1121 
show the predicted and observed L/E distributions for 
these two samples, before and after the separation of 
events into neutrinos and antineutrinos. For each of these 
figures, the observed data are compared with the best fit 
neutrino and antineutrino oscillation parameters, which 
are given in Section [Villi 
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FIG. 9. Distributions of track fit uncertainty, {q/p)/cr q / p and track straightness variable, Xune/dof, used to select events with 
well- measured muon charge sign. The distributions are plotted for contained-vertex muons (top panels), and neutrino-induced 
rock-muons (bottom panels). In each plot, the dashed line indicates the total prediction in the absence of oscillations; the 
solid line shows the prediction for oscillations with Am 2 = 2.32 x 10~ 3 eV 2 and sin 2 26 l = 1.0; the shaded histogram shows the 
cosmic-ray muon background; and the points show the observed data. In addition, the hatched histograms show the component 
with mis-identified charge sign. The arrows indicate the selections used to identify events with well-measured charge sign. 
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FIG. 10. Distributions of reconstructed zenith angle, for contained-vertex muons (left), and neutrino- induced rock-muons 
(right). In each plot, the dashed line gives the nominal prediction for the case of no oscillations; the shaded histogram shows 
the cosmic-ray muon background; and the points with errors show the observed data. The solid line shows the best fit to the 
data, which combines the best fit oscillation and systematic parameters. 
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FIG. 11. Distributions of reconstructed log 10 (L/_E), for contained-vertex muons (left), and neutrino-induced rock-muons (right). 
In each plot, the dashed line gives the nominal prediction for the case of no oscillations; the shaded histogram shows the cosmic- 
ray muon background; and the points with errors show the observed data. The solid line shows the best fit to the data, which 
combines the best fit oscillation and systematic parameters. 



Contained-vertex muons Qi-) 




co 

-»— * 
> 

LU 



40 - 



20 - 



log io (L[km]/E[GeV]) 
Contained-vertex muons 



V 



I 1 1 1 1 I 1 1 — 1 ' I 1 ' — ' 1 I 

No Oscillations 
MINOS Best Fit 
Cosmic-ray jx 

Data 




12 3 4 
log io (L[km]/E[GeV]) 



Neutrino-induced rock-muons pi-) 




30 h — 

V 



£ 20 - 
c= 

> 
UJ 

10 



log io (L[km]/EJGeV]) 
Neutrino-induced rock-muons 



~i ■ — ■ ■ ■ i ■ ■ ■ 1 i 
No Oscillations 
MINOS Best Fit 



• Data 




log io (L[km]/EJGeV]) 



FIG. 12. Distributions of reconstructed log 10 (L/E), plotted for selected neutrinos and antineutrinos in the contained-vertex 
muon sample (top panels) and neutrino-induced rock- muon sample (bottom panels). In each case, the dashed line gives the the 
nominal prediction in the absence of oscillations; the shaded histogram shows the cosmic-ray muon background; and the points 
with errors show the observed data. The solid line indicates the best fit to the data, which combines the best fit oscillation and 
systematic parameters. 
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A. Separation of events by L/E resolution 

The intrinsic L/E resolution of contained- vertex 
events, and therefore the degree to which they contribute 
to the overall sensitivity to the oscillation parameters, 
varies significantly across the high resolution sample. 
The resolution in the propagation distance depends on 
the reconstructed energy and zenith angle, and is worse 
at low energies where the average angle between the neu- 
trino and muon is large, and also around the horizon 
where the propagation distance varies rapidly as a func- 
tion of zenith angle. The resolution in the neutrino en- 
ergy is worse for events where the muon momentum is 
determined from curvature rather than range, and also 
for events with high y values, since the shower energy 
resolution is generally poorer than the muon momentum 
resolution. 

The sensitivity to oscillations can be improved by in- 
corporating information on L/E resolution into the oscil- 
lation fit. For this analysis, a Bayesian technique is used 
to estimate the L/E resolution of selected contained- 
vertex muons on an event-by-event basis. For each event, 
a probability distribution function (PDF) in \og 10 (L/E) 
is calculated by combining the measured muon momen- 
tum, muon direction and shower energy of the event with 
information from the Monte Carlo simulation describ- 
ing the atmospheric neutrino spectrum, interaction kine- 
matics and detector resolution. The simulation is used 
to construct PDFs relating the measured muon momen- 
tum and shower energy of selected i/ M and F M events to 
their corresponding true distributions. The simulation 
also provides PDFs in the kinematic variables W 2 and y 
for and v ^ CC interactions, binned as a function of 
neutrino energy, which enable the distributions of muon 
momentum and shower energy to be mapped onto a dis- 
tribution of neutrino energy, and the muon direction to 
be mapped onto a distribution of neutrino zenith angle. 
The overall PDF in \og 10 (L / E) is then obtained by tak- 
ing a convolution of these neutrino distributions and the 
RMS of this PDF, (Ti og (L/E), gives the L/E resolution. 
A full description of the technique is given in |47| . 

Figure [13] shows the predicted and observed <Ji og (L/E) 
distributions for high resolution contained-vertex muon 
neutrinos. The shape of the predicted distribution is al- 
most independent of the input oscillation parameters. 
The observed spread of <Ji og (L/E) values is substantial, 
and corresponds to 25% of the spread in \og w {L/E). 
Therefore, a significant gain in sensitivity is expected by 
separating events into bins of L/E resolution. 

For the oscillation analysis, the selected contained- 
vertex muon neutrinos are divided into the following four 
bins of L/E resolution: 

0.00 < a log(L/E) < 0.25, 
0.25 < a log(L/E ) < 0.50, 
0.50 < <J log(L/E) < 0.75, 
0-75 < a log{L/E) < 1.50. 

Figure Q3] shows the ratio of the predicted L/E distri- 



butions with oscillations to those without oscillations in 
each bin of resolution. The oscillations are most sharply 
resolved in the bin with the best L/E resolution. Here, 
the ratio initially falls with L/E, reaching a minimum at 
the peak oscillation probability. The ratio subsequently 
rises to a maximum, and a second oscillation dip is visible 
before the ratio averages to 1 — isin 2 2# as the frequency 
of oscillations becomes rapid. 

Selected neutrino-induced rock-muons are separated 
into low momentum {P^ < 10 GeV) and high momen- 
tum (P M > 10 GeV) samples. This separation roughly 
distinguishes those muons whose parent neutrinos have a 
relatively large oscillation probability from those with a 
lower probability. 

Figure [T5l shows the predicted and observed L/E dis- 
tributions, separated into bins of L/E resolution for 
contained-vertex muons, and into bins of muon momen- 
tum for neutrino-induced rock-muons. The predicted dis- 
tributions are calculated for the case of no oscillations, 
and for the best fit neutrino and antineutrino oscillation 
parameters. 

B. Oscillation fit 

An oscillation fit is applied to the data assuming two- 
flavor — > v T vacuum oscillations. In this approxima- 
tion, the oscillation probability is given by: 

9 n 9 /Am 2 L\ 
P -> v T ) = sin 2 20 sin 2 ( \ , 

where Am 2 and sin 2 29 are the two- flavor oscillation pa- 
rameters, L is the neutrino propagation distance, and E 
is the neutrino energy. 

For upward-going atmospheric neutrinos with energies 
in the 2 — 10 GeV region, an asymmetry between muon 
neutrinos and antineutrinos is predicted, arising from the 
resonant enhancement of three-flavor oscillations by mat- 
ter effects [48|. The sign of the asymmetry depends on 
the sign of the neutrino squared-mass difference, and is 
therefore sensitive to the mass ordering of neutrinos [49| . 
These effects have only a small influence on the predicted 
L/E distributions of neutrinos and antineutrinos and so 
are not considered in this analysis. 

The high and low resolution contained-vertex muons, 
neutrino-induced rock-muons, and contained-vertex 
showers are each included as separate samples in the 
oscillation fit. For the low resolution contained-vertex 
muons and contained-vertex shower samples, the events 
are fitted in single bins of normalization. The high reso- 
lution contained-vertex muon sample is divided into two 
bins of direction, r — (u,d), corresponding to upward- 
going (u) and downward-going (d) muons; three bins of 
charge sign, s=(u,V,X), corresponding to neutrinos (u), 
antineutrinos (V) and events with ambiguous charge sign 
(A); and four bins of L/E resolution. The neutrino- 
induced rock-muon sample is divided into three bins of 
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FIG. 13. Distributions showing the calculated L/E resolu- 
tion for the high resolution sample of contained-vertex muon, 
which have well-measured propagation direction. The dashed 
histogram indicates the nominal prediction in the absence of 
oscillations; the solid histogram indicates the prediction for 
the best fit oscillation parameters presented in this paper; 
the shaded histogram shows the cosmic-ray muon background; 
and the points with error bars represent the observed data. 
The vertical dashed lines correspond to the partitions used to 
divide the selected events into four bins of L/E resolution. 
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FIG. 14. Ratios of the predicted \og 10 (L/E) distributions 
with oscillations to those without oscillations in the four bins 
of L/E resolution. The predictions with oscillations are gen- 
erated using input parameters of Am 2 = 2.32 x 10 -3 eV 2 and 
sin 2 28 = 1.0. The oscillations are most sharply defined in 
the bin of highest resolution. Here, a clear oscillation dip can 
be seen at \og 10 (L/E) ~ 2.7, corresponding to the peak os- 
cillation probability. The ratio then rises to a maximum at 
\og 10 (L / E) w 3, and a second dip is visible before the ratio 
averages to 1 — isin 2 2# = 0.5, as the frequency of oscillations 
becomes rapid. 



charge sign, and two bins of muon momentum. Overall, 
there are 24 high resolution contained-vertex muon and 
6 neutrino-induced rock-muon distributions. Each dis- 
tribution is binned in log w (L/E), using 25 bins in the 
range [—0.5, +4.5]. This gives a total of 750 high resolu- 
tion bins, in addition to the 2 low resolution bins. 

A maximum likelihood fit to the data is performed us- 
ing the following negative log-likelihood function: 



In C — fx — n In li + fi — n In /i 

l h r,s 



h r,s i.k 
v 2 



E 



07 



2al. 

This log-likelihood function is divided into the following 
terms: 

1. Normalization: The sums ^2 f- — n m M represent 
the Poisson probability for observing a total of n 
events with a prediction of \x events. The first sum, 
denoted I, is taken over the contained-vertex shower 
and low resolution contained-vertex muon samples, 
which are fitted in single bins; the second sum, de- 
noted h, is taken over the neutrino-induced rock- 
muon and high resolution contained-vertex muon 
samples, which are separated by muon direction 
r = (u, d) and charge sign s = (y, F, X). 

2. Shape Term: The shape of the log 1Q (L/E) distri- 
bution is incorporated into the oscillation fit for 
the neutrino-induced rock-muon and high resolu- 
tion contained-vertex muon samples. The terms 
Yli k Uik m (fik) represent the likelihood functions 
for each of the \og 10 (L/E) distributions included 
in the fit. The i-sum is taken over each resolu- 
tion bin for the contained-vertex muons, and each 
momentum bin for neutrino-induced rock-muons; 
the fc-sum is taken over each of the 25 bins in the 
log 10 (L/E) distribution. Within the sum, is 
the observed number of events and fik is the rela- 
tive predicted probability in the i th and k th bins. 

3. Systematic Uncertainties: Systematic effects are 
incorporated as nuisance parameters, where the 
shift ctj is the deviation of the j th systematic pa- 
rameter from its nominal value. A penalty term, 
ev 2 /2<7^, , is added to the likelihood, where the er- 
ror a aj represents the estimated uncertainty in the 



■th 



systematic parameter. 



A total of 12 systematic uncertainties are incorporated 
into the fit as nuisance parameters, as listed in Table Hill 
For contained-vertex neutrino interactions, a 15% uncer- 
tainty is applied to the normalization of the event sam- 
ple. The following additional uncertainties are applied 
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FIG. 15. Distributions of log 10 (L/ E) observed in each bin of L/E resolution for contained-vertex muons, and in each bin of 
muon momentum for neutrino-induced rock-muons. For each of the panels, the dashed line gives the the nominal prediction in 
the absence of oscillations; the shaded histogram shows the cosmic-ray muon background; and the points represent the data. 
The solid line indicates the best fit to the data, combining the best fit oscillation and systematic parameters. 



to this sample: a 3% uncertainty on the up-down ratio; 
a 5% uncertainty on the {v^ 1 +T' fl )/{v e +T> e ) ratio; a 10% 
uncertainty on the v^jv^ ratio; and a 20% uncertainty 
on the ratio of NC to CC interactions. For neutrino- 
induced rock-muons, a 25% uncertainty is applied on the 
normalization of the event sample. An additional 12.5% 
uncertainty is applied to the v ^jv^ ratio in this sample. 

To account for the uncertainty in the shape of the 
atmospheric neutrino energy spectrum, the number of 
contained-vertex events is allowed to scale as a function 
of neutrino energy. The form of the scaling function is 
chosen to cover the variations in the spectrum generated 
by changing the flux model and by reweighting the cross- 
section model according to its given uncertainties. Above 
3 GeV, where the prediction from the simulation approx- 
imately follows a power function, events are scaled by 
f{E v ) = 1 + a \n{E v /?>). Below 3 GeV, this is connected 
smoothly to a linear function, f{E u ) = 1 + a (E v — 3). 
The scaling function is applied separately to neutrinos 
and antineutrinos. In each case, the spectrum parameter, 
a, is normally distributed with a standard deviation of 
6%. Finally, to account for the systematic uncertainties 
on the track and shower energy scale, a 3% uncertainty 
is included on the muon momentum from range, 5% on 
the momentum from curvature, and 15% on the shower 
energy scale. Of all the systematic uncertainties incorpo- 



rated into the fit, only the two normalization parameters 
are found to have a significant impact on the resulting 
confidence limits. 



C. Results of oscillation fit 

The log-likelihood function is minimized with respect 
to the oscillation and nuisance parameters. Table IIIII 
summarizes the best fit parameters. The best fit point 
occurs at (|Am 2 |, sin 2 26») = (1.9 x 10~ 3 eV 2 , 0.99). The 
68%, 90% and 99% confidence limits (C.L.) on the oscilla- 
tion parameters are obtained in the limit of Gaussian er- 
rors from the locus of points with log-likelihood values of 
-Aln£ = (1.15, 2.30, 4.61) relative to the best fit point. 
Figure [TB] shows the resulting 90% contours from this 
analysis. For comparison, this figure also shows the 90% 
contours from the MINOS beam neutrino analysis [Hj], 
and also from the Super-Kamiokande atmospheric neu- 
trino zenith angle analysis (2lj . 

The log-likelihood surface is used to calculate single- 
parameter confidence intervals for each of the oscilla- 
tion parameters, by minimizing with respect to the other 
oscillation parameter. The 90% single-parameter con- 
fidence intervals at the best fit point, calculated using 
this method, are |Am 2 | = (1.9 ± 0.4) x 10~ 3 eV 2 and 
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Parameter 


Uncertainty 


Best Fit (2 Osc. Param.) 


Best Fit (4 Osc. Param.) 


|Am 2 |/eV 2 
|Am 2 |/cV 2 
sin 2 26» 






1.9 x 10~ 3 


2.2 x 10~ a 






1.9 x 10~ 3 


1.6 x 10~ 3 






0.99 


0.99 


sin 2 26 






0.99 


1.00 


normalization (contained-vertex v) 


a - 


= 15% 


+0.6 a 


+0.7a 


normalization (^-induced rock-pi) 


a ■ 


= 25% 


+0.1 a 


+0.1 a 


up /down ratio (contained-vertex v) 


a 


= 3% 


-0.1a 


-0.1a 


Vejv^ ratio (contained-vertex v) 


a 


= 5% 


-0.5 a 


-0.5 a 


v^/Vfj, ratio (contained-vertex v) 


a - 


= 10% 


-0.5 a 


-0.6a 


v^jvp ratio (f -induced rock-/i) 


a = 


= 12.5% 


+1.1 <7 


+0.9 a 


NC/CC ratio (contained-vertex v) 


a ■ 


= 20% 


+0.6a 


+0.6a 


v spectrum parameter 


a 


= 6% 


-0.4a 


-0.4a 


V spectrum parameter 


a 


= 6% 


+0.3 a 


+0.3 a 


u momentum (range) 


a 


= 3% 


-0.3a 


-0.3 a 


jx momentum (curvature) 


a 


= 5% 


+0.3 a 


+0.3 a 


shower energy 


a - 


= 15% 


+0.4 a 


+0.4 a 



TABLE III. Summary of systematic uncertainties included in the oscillation fit, along with the best fit oscillation and systematic 
parameters returned by each fit. For the two-parameter fit, equal oscillation parameters are used for neutrinos and antineutrinos; 
for the four-parameter fit, separate oscillation parameters are used. The best fit systematic parameters are given in units of 
standard deviations. 
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FIG. 16. Confidence limits on the parameters |Am 2 | and 
sin 2 2f9, assuming equal oscillations for neutrinos and antineu- 
trinos. The solid line gives the 90% contour obtained from 
this analysis, with the best fit parameters indicated by the 
star. For comparison, the dashed line shows the 90% contour 
given by the MINOS oscillation analysis of neutrinos from 
the NuMI beam [Q, with the best fit point indicated by the 
triangle. The dotted line shows the 90% contour from the 
Super-Kamiokande atmospheric neutrino zenith angle analy- 
sis (from [10]]), with the best fit point indicated by the circle. 



sin 29 > 0.86. The null oscillation hypothesis is disfa- 
vored at the level of 9.2 standard deviations. 



VIII. FITS TO NEUTRINO AND 
ANTINEUTRINO OSCILLATIONS 

Since the data are separated into pure samples of neu- 
trinos and antineutrinos, they can be used to study os- 
cillations separately in neutrinos and antineutrinos. The 
analysis described above is extended to incorporate sep- 
arate oscillation parameters for neutrinos (Am 2 , sin 2 20) 
and for antineutrinos (Am 2 , sin 2 20). The log-likelihood 
function is then minimized with respect to these oscilla- 
tion parameters and the twelve nuisance parameters. The 
best fit occurs at (|Am 21 " : " J " 

„ 5 „2r 



1,81^20) = (2.2 X 10- 3 eV^,0.99) 
and (|Am 2 |, sin 2 20) = (1.6 x 10~ 3 eV 2 , 1.00), as given 
in Table IIIII The neutrino and antineutrino oscillation 
parameters are found to be approximately uncorrelated 
around the best fit point. A set of two-parameter profiles 
can be calculated from the four-parameter likelihood sur- 
face by minimizing with respect to pairs of oscillation pa- 
rameters. Figure IT71 shows the resulting 90% contours ob- 
tained for the (|Am 2 |, sin 2 20) and (|Am 2 |, sin 2 20) planes. 
These results are compared with the 90% contours from 
the MINOS analyses of NuMI beam data acquired in neu- 
trino [12J and antineutrino [15J mode, and also the 90% 



contours from the SK analysis of atmospheric neutrinos 
and antineutrinos (2lj . 

The four-parameter likelihood surface is used to calcu- 
late single-parameter confidence intervals on each of the 
four oscillation parameters. The resulting 90% C.L. are: 
lO" 3 eV 2 and sin 2 20 > 0.83 for neu- 

r2 



| Am 2 1 = 2.2t 2 A 6 x 

trinos; and |Am 2 | = 1.6±gj x 10~ 3 eV 2 and sin 2 > 0.76 
for antineutrinos. The null oscillation hypothesis is disfa- 
vored at the level of 7.8 standard deviations for neutrinos 
and 5.4 standard deviations for antineutrinos. 

As a measure of the quality of the fit, a set of 10,000 
simulated experiments were generated at the best fit os- 
cillation parameters. For each simulated experiment, in- 
put systematic parameters were chosen from Gaussian 
PDFs with widths set to the systematic uncertainties. 
The best fit parameters were then found for each experi- 
ment by minimizing the log-likelihood function. For each 
experiment, the minimum value of — A In C was recorded; 
in 22% of experiments, the value exceeded that obtained 
from the fit to the data. 

Figure [LSI compares the observed 90% C.L. from each 
fit with the predictions from the Monte Carlo simulation, 
calculated by inputting the best fit oscillation parameters 
into the simulation. For the two-parameter oscillation fit, 
where neutrinos and antineutrinos take the same oscil- 
lation parameters, there is good agreement between the 
observed and predicted contours. For the four-parameter 
oscillation fit, where neutrinos and antineutrinos take 
separate oscillation parameters, there is a good match 
between contours for the limits on the sin 2 29 and sin 2 29 
parameters and the lower limits on the |Am 2 | and |Am 2 | 
parameters. However, the upper limits on these parame- 
ters are found to be higher than predicted for neutrinos 
and lower than predicted for antineutrinos. 

As a check on the observed confidence limits, the full 
likelihood surface was calculated for a set of 250 simu- 
lated experiments, generated at the best fit oscillation 
parameters from the two-parameter fit. The resulting 
90% confidence intervals were then calculated for each 
experiment. In 25% of these experiments, the confidence 
intervals obtained for the Am 2 parameter are broader for 
neutrinos than antineutrinos, as is the case for the ob- 
served data; in 10% of the experiments, the relative size 
of these intervals is larger than for the observed data. 
These results indicate that the confidence intervals cal- 
culated from the observed data are reasonable. 

Finally, a log-likelihood profile is calculated in the 
(| Am 2 1, | Am 2 1) plane, by minimizing the log- likelihood 
function with respect to the sin 2 20 and sin 2 20 param- 
eters. Figure Qj5] shows the resulting 68%, 90% and 
99% confidence intervals. This log-likelihood profile is 
used to place limits on the difference between the neu- 
trino and antineutrino mass splittings |Am 2 | and |Am |. 
The single-parameter 90% confidence intervals, assuming 
Gaussian errors, are |Am 2 | — |Am 2 | = 0.6+g'g x ^ 3 e V 2 - 
This result is consistent with equal mass splittings for 
neutrinos and antineutrinos. 
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FIG. 17. Confidence limits on neutrino (left) and antineutrino (right) oscillation parameters. The solid lines show the 90% 
contours obtained from this analysis. The two-parameter contours for neutrinos and antineutrinos are calculated by profiling 
the four-parameter likelihood surface. The best fit parameters are indicated by the stars. For comparison, the dashed lines in 
each plot show the 90% contours from the MINOS analysis of beam data in neutrino mode [12| and antineutrino mode [15|, 
with the best fit points indicated by the triangles. The dotted lines show the 90% contours from the Super-Kamiokande analysis 
of atmospheric neutrinos and antineutrinos (from |2lj]), with the best fit points indicated by the circles. 




FIG. 18. Comparisons of the observed and predicted 90% contours. The left plot shows the contours for the two-parameter 
oscillation fit, where neutrinos and antineutrinos take the same oscillation parameters. The right two plots, labelled v and 
V, show the contours obtained for neutrinos and antineutrinos, respectively, resulting from the four-parameter oscillation fit, 
where neutrinos and antineutrinos take different oscillation parameters. In each case, the predicted contours are generated by 
inputting the best fit parameters from the observed data into the simulation, and running the full oscillation fit, including the 
twelve systematic parameters. 
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FIG. 19. Confidence limits obtained for the oscillation pa- 
rameters | Am 2 1 and |Am 2 j, representing the mass splittings 
for neutrinos and antineutrinos, respectively. At each point 
in parameter space, the negative log-likelihood function has 
been minimized with respect to the mixing parameters sin 2 29 
and sin 2 26. The 68%, 90% and 99% contours are indicated 
by the dotted, solid and dashed curves, respectively, and the 
best fit parameters are indicated by the star. The diagonal 
dashed line indicates the line of IAm 2 | = |Am 2 |. 



spheric neutrino events. The events are separated into 
905 contained-vertex muons and 466 neutrino-induced 
rock-muons, produced by and F M CC interactions, 
and 701 contained-vertex showers, composed primarily 
of v e and v e CC interactions and NC interactions. The 
curvature of muons in the magnetic field is used to di- 
vide the selected contained-vertex muons and neutrino- 
induced rock-muons into separate samples of and V y, 
events. The double V^jv^ ratio is calculated to be: 
^ata/RMC = 1 03 ± q 08 (stat.) ± 0.08 (syst.). 

A maximum likelihood fit to the observed L/E distri- 
butions is used to determine the atmospheric neutrino 
oscillation parameters. The sensitivity to oscillations is 
improved by separating the contained-vertex muons into 
bins of L/E resolution, and the neutrino- induced rock- 
muons into bins of muon momentum. The fit returns 



90% confidence limits of |Ar 



(1.9 ±0.4) x 10-^eV 



and sin 2 2# > 0.86. The oscillation fit is extended to 
allow separate oscillation parameters for neutrinos and 
antineutrinos. This fit returns 90% confidence limits of 
|Am 2 | - |Am 2 | = O.eto.s x 10~ 3 eV 2 on the difference 
between the squared-mass splittings for neutrinos and 
antineutrinos, consistent with equal mass splittings for 
neutrinos and antineutrinos. 
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